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ABSTRACT 

We conduct a systematic study of the angular momentum problem in numerical simu- 
lations of disk galaxy formation. We investigate the role of numerical resolution using 
a semi-cosmological setup which combines an efficient use of the number of parti- 
cles in an isolated halo while preserving the hierarchical build-up of the disk through 
the merging of clumps. We perform the same simulation varying the resolution over 4 
orders of magnitude. Independent on the level of resolution, the loss of angular momen- 
tum stays the same and can be tied to dynamical friction during the build-up phase. 
This is confirmed in a cosmological simulation. We also perform simulations including 
star formation and star formation and supernova feedback. While the former has no 
influence on the angular momentum problem, the latter reduces the loss to a level 
potentially in agreement with observations. This is achieved through a suppression 
of early star formation and therefore the formation of a large, slowly rotating bulge. 
We conclude that feedback is a critical component to achieve realistic disk galaxies in 
cosmological simulations. Numerical resolution is important, but by itself not capable 
of solving the angular momentum problem. 

Key words: galaxies: spiral - formation - evolution - methods: iV-body simulations 
- hydrodynamics 



1 INTRODUCTION 



' The prevailing picture of gal axy formation was fo rmu- 
. lated almost 30 years ago by IWhite fc ReesI l| 19781) and 
. iFall fc Efstathioul (Il980l ). Gas falls into the potential wells 
' of hierarchically growing dark matter structures and is ad- 
ditionally governed by dissipational processes. The idea of 
dark matter halos acquiring angular momentum through 
tidal torques fro m surrounding lar ge sc ale structure was 
first developed by jStromberd lll933) and iHovld lll94£ ) and 
later quantified by Peeblef^ ('1969') , 'Doroshkevich| (|l97oi ) and 
[white (1984). It results in slowly rotating dar k matter halo s 
iBarnes fc Efstathioull 19871 : ISteinmetz fc Bartelmannlll995h . 



The gas within the halo acquires the same specific angu- 
lar momentum, but it can cool, therefore collapses and, 
under of cons ervation of angular momentum, forms disks 
||Mo. Mao fc W hite 1998). Stars are formed due to the frag- 
mentation of the gas disks, and spheroi ds are formed dur- 
ing th e interaction and merging of disks (jToomre fc Toomrd 
1 19721 ). The detailed interaction of gas, stars and related 
physical processes make this a very complex problem, which 
needs to be addressed through numerical simulations. This 
is a twofold challenge. The first challenge is posed by the 
large dynamic range of scales which are encountered, from 
the tens of Mpc scales of halo and galaxy interaction during 
assembly to the sub-pc scales of star formation. The second 



challenge stems from physical processes like star formation 
and feedback, which happen on scales below the resolution 
of most present day cosmological simulations. It is thus not 
so surprising that, despite many efforts in the past years, 
disk galaxies with realistic properties remain elusive in simu- 
lated universes. Especially bulge-less disks, w hich have been 
show n to exist in non-negligible numbers l|Kautsch et al.l 
[2006') in the nearby universe, have not been produced in 
a simulation up to this point. 

One of the major problems encountered in nu- 
merical simulations of disk for mation is the so calle d 
angular momentum catastrophe (|Navarro fc Benj Il99ll ). 
Early simulations all resulted in small, highly centrally 
conc entrated disks w i th dominating bulges (for exam- 
ple iNavarro fc White! Il994l : ISteinmetz fc Navarro! !l999! : 
!Navarro fc Steinmet j !2000! ) after loosing large amounts 
of angular momentum during the formation. As four 
main causes have been proposed: (1) lack o f resolution 
ijGovernato et al.!!200j : !Kaufmann et al.!!2007l ): (2) an ar- 
tificial transfer of angular mom entum from the cold disk 
at the interface to the hot halo ( jOkamoto et al."200 5j); (3) 
disk instabilities leading to fragmentation (jRobcrts on et al.! 
!2004 ) ; (4) merging of substructure leading to transfer of an- 
gular momentum to the outer halo via dynam ical friction 
(iNavarro fc Benj!l99l! : !Navarro fc White!!l994! ) . The second 
and third effect can be avoided with an improved imple- 
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mentation of gas physics (a decoupling of h ot and cold gas 
hase s for the second effect as suggested bv lOkamoto et ahl 
or a multiph ase medium for the third effect according 



Robertson et al. I ([2004 )). Low numerical resolution, the 



to 

first effect, can cause artificial loss of angular momentum 
through a variety of effects like numerical two-body heating 
due to the la rge mass difference betwe e n dark matt er and 
gas particles (|Steinmetz fc White! [l997l : iMaverl 120051 ). arti- 
ficial viscosity, an incorrect pressure estimate in the stan- 
dard SPH method for particles with significant tempera- 
ture differences (as happens for example at the edge of the 
disk) , and an artificial asymmetry of the disk causing gravi- 
tational torques from the surrounding halo. These problems 
can be remedied with larger numbers of particles in the disk. 
Others, like a correct modeling of the disk center where ve- 
locity gradients are steep and angular momentum is easily 
lost, or the formation of bars, which can transport angular 
momentum outward efficiently, are difficult to address even 
with very large particle numbers. While these effects have 
been demonstrated to oper ate in well established d isks and 
smooth collapse scenarios (|Kaufmann et all |2007^ . it still 
needs to be demonstrated to what extend they are responsi- 
ble for the angular momentum loss seen in fully cosmological 
large scale structure simulations. 

As opposed to the first three effects, the fourth, merg- 
ing of substructure, is a physical effect which is natural to 
ACDM. ft affects dark matter as well as gas, and angular 
momentum is lost through dynamical friction of interact- 
ing and merging clumps. The most efficient mechanism to 
prevent angular momentum loss through this process is to 
prevent the formation of too many clumps. For this, two sce- 
narios have been discussed. The first one involves a change 
in the cosmological model from cold to warm dark matter 
(WDM). While leaving the large scales untouched, this re- 
sults in a smoothing of small scales of the power spectrum, 
and therefore less substructure and a smoother accretion 
of baryons onto dark matter hal os instead of dense clumps 
suscep tible to dynamical friction. ISommer-Larsen fc Dolgovl 
([2OOII) successfully showed that this leads to an alleviation 
of the angular momentum problem, which puts doubt on 
the claim that the angular momentum problem is mainly 
a resolution issue. It can also help to solve what is known 
as the "missing satellite problem" in ACDM, a large excess 
in dark matter substructure predicted by numerical simu- 
lations compared to observed numbers of satellite galaxies 
l|Gotz fc Sommer-Larsenll2003l ). However, observational con- 
straints on the mass of dark matter particles, for example 
by the redshift of re-ionization (which is delayed significantly 
in a WDM universe due to the lack of low mass halos), al- 
ready allow only varieties which are tepid at most l|Viel et al.l 
|2008| ). In the second scenario preventing the gas from form- 
ing small, dense clumps before entering the main halo, the 
thermal energy of the gas is increased via physical processes 
like feedback from supernova explosions or stellar winds. We 
will address this in the second part of this paper. 

In this paper, we will focus on the role of resolution and 
feedback with respect to the angular momentum problem in 
numerical simulations of disk galaxy formation. We present 
a detailed, systematic study of resolution effects in isolated 
but semi-cosmological halos preserving the essence of the 
hierarchical formation history, and also test the results in 
fully cosmological runs. The study also covers the effects of 



various physical processes ranging from simulations includ- 
ing only radiative cooling over the inclusion of simple star 
formation to the application of a model for supernova feed- 
back. The paper is organized as follows: in §2 we describe the 
code and the initial conditions we use. In §3 we present the 
results on the angular momentum problem in simulations in- 
cluding only gravity and radiative cooling. §4 discusses the 
impact of star formation and the feedback model and in §5 
we summarize and discuss our results. 



2 CODE AND INITIAL CONDITIONS 

All simulations ar e based on the WMAP3 cosmology 
ISpergel et al.ll2007D with Ho=73 km s-^Mpc-\ as = 0.75, 
ns = 0.9, ^0 = 0.24, = 0.76 and »b = -04. 
We use the N-body code GADGET2 jSpringell l2005h in 
a version including radiative cooling based on iKatz et al.l 
(1996). To minimize the number of free parameters we do 
not include a UV background. The gas is treated within 
the Smoothed Particle Hydr o dynamics (S PH) framework 
dGingold fc MonaghanI l| 19771 ). ILucvI (|l977l ')). In the follow- 
ing, we explain our implementation of star formation and 
supernova feedback. 



2.1 Star formation and feedback 

Our basic star fo r mation model follows iKat j l| 19921 ) and 
iNavarro fc Whitel (|1994| ). The main criteria for either 
spawning a new star particle from a gas particle or turn- 
ing the gas particle into a star particle (depending on the 
gas mass) are: 

(1) p > pcrit, where we chooseQ /9crit=7x 10~^^g cm"'' 

(2) < 

(3) T< 3 X 10* K. 

The star formation rate is guided by a basic Schmidt law 
(Schmidt 1959.), but expressed in physical density as 

dp^^ (1) 
dt U 

where c* is a dimensionless parameter regulating the star 
formation effi ciency (we u s e c<-= 0.1, in agreement with ob- 
servations by iDuerr et all (jl982h for the ISM in the Milky 
Way) and 

t^, = max{tdyn,tcool), tdyn = (47rGp) ~ ^''^ . (2) 

We use a stochastic approach to star formation, so the prob- 
ability to form a star particle during a given timestep At is 
given by 



P = 



(3) 



e is a parameter giving the fraction of the gas mass turned 
into a star or, alternatively interpreted, the number of gen- 
erations of stars. We generally use e—2 to prevent too much 
increase in the number of particles. 



^ This choice of density corresponds to a region where the cool- 
ing time is always shorter than the dynamical time, resulting 
in rapid cooling to the cut-ofE temperature of the cooling curve 
(lO^'K) where the Jeans instability criterion (sound crossing time 
is larger than the dyn amical time) is almost always fulfilled 
iNavarro fc WhitdlToM) . 



The Angular Momentum Problem in Cosmological Simulations of Disk Galaxy Formation 3 



We include the most basic but also probably most influ- 
ential type of stellar feedback: energy output from the late 
stages of the evolution of massive stars (M>8 M©). This 
is usually referred to as supernova feedback, but actually 
covers all types of energetic feedback which are not resolved 
in this study. The model implemented here showed the best 
performance in a m ore extensive study of feedba ck models 
which is presented in lPiontek fc Steinmet j (|2009bl ). As in all 
N-body simulations, each "star particle" represents a whole 
population of stars of total mass (which depends on res- 
olution). For this we a ssume a Miller-Sc alo initial mass func- 
tion (IMF) (IMF, Mi ller fc Scald 1 19791 ) with a lower cutoff 
of 0.1 M0 and an upper cutoff of 100 M©: 



C(M) = M^A 




0.1 < M < IM0 

1 < M < 2 M0 

2 < M < IOM0 
10 < Af < 100 Mq 



(4) 



where A=0. 284350751 and M* is the mass of the stellar 
population, in our case of the star particle. Stars with masses 
between 8 and 100 Mq explode as type II supernova and 
each explosion is assumed to yield an energy of 10^^ ergs. 
This corresponds to an energy of 1.21x10^® erg per solar 
mass formed. To avoid numerical artefacts related to abrupt 
changes in the thermal properties of the gas which cannot 
be resolved, the energy is returned to the surrounding gas 
over time following an exponential approach using 

AEit) = EsN—^e~~—. (5) 

tsN tsN 

Here, Esn is the total supernova energy, t* is the time when 
the star particle was created and tsN is a typical time scale 
(we use 20 Myr as the lifetime of a 8 Mq star). This energy is 
smoothed over the neighboring gas particles using the SPH 
smoothing kernel, so each neighbor gets an energy of 



AE{t)- 



(6) 



h* is the smoothing length and the most important parame- 
ter here. This is a variable length scale calculated similar to 
the gas smoothing length in GADGET2 (the mass within 

the smoothing s pher e is kept constant). 

iKatj l|l992l 'l and lSteinmetz fc Miilleil (liggsl ) found that 
supernova feedback energy is radiated away very quickly in 
the simulations, since star formation by definition occurs 
in dense regions with short cooling times. Consequently, 
in these early simulations, feedback had little impact. In 
more realistic ISM models this does not happen due to the 
multiphase nature of the gas, but it is computationally de- 
manding, if at all possible, to achieve enough resolution to 
directly simulate this. One meth od to prev e nt thi s immedi- 
ate re-radiation was suggeste d bv lGerritsenI (Il997l ) and later 
iThacker fc Couchmanll2000l ). It involves turning off radia- 
tive cooling in gas particles receiving feedback energy and 
evolving them adiabatically for a given time. We adopt this 
method and turn off the cooling locally for a time of 20 Myr 
(in agreement with earlier simulations and similar to the 
lifetime of a supernova shell). 

2.2 Isolated initial conditions 



simulation, this is difficult, since the resolution for a given 
halo and particularly for the forming galaxy is limited and 
the hierarchical merging processes are likely to destroy disk 
progenitors. Therefore we first perform a controlled exper- 
iment for an isolated halo, and then verify our conclusions 
in a cosmological box. However, a completely isolated and 
idealized setup consisting o f an NFW halcH in combination 
with a hot gas halo, as in iKaufmann et al.l (|2007l '). is too 
artificial to probe the relevant physics, since it completely 
lacks any merging in the formation process. To be more re- 
alistic with the same benefits of high resolution in the gas 
disk and no major mer gers destroy ing it, we use a semi- 
cosmological setup as in iKat j (|l99lh . We start with a cos- 
mological box with a side length of 2 h~^Mpc, large enough 
to host one Milky Way-type halo, but lacking all large scale 
structure. Centering the forming halo, we cut out a sphere 
with 1 h^^Mpc radius which is simulated in physical coordi- 
nates and with vacuum boundary conditions. To mimic the 
missing large scale effects, we apply an overdensity of 3a 
and put the sphere into solid body rotation with a spin pa- 
rameter of A = 0.08. This is high er than what is typical for 
dark matter halos (A f» 0.04 B ullock et al.ll200lf ). but will 
ensure the formation of a large rotationally supported and 
stable gas disk. We include gas particles by splitting each 
dark matter particle into a gas and a dark matter particle 
with the mass ratio given by the density parameters. The 
initial redshift of the simulations is z=75. For our resolu- 
tion study, we create the same initial conditions with four 
different resolutions, based on boxes with 16^, 32^, 64^ and 
128"^ dark matter particles. This corresponds to spheres with 
2 X 10^, 2 X 10'*, 1.4 X 10^ and 1.1 X lO'' particles each in gas 
and dark matter. We use a gravitational softening growing 
proportional to the scale factor for z>10 and staying fixed 
in physical coordinates after z=10, at a value of lh~^kpc 
for the 3 lower resolutions and 0.25 h~'^kpc for the highest 
resolution, for both gas and dark matter particles. We also 
tested higher softenings for the two lowest resolutions - 4 
kpc h"^ for 16^ and 2 kpc h'^ for 32 ^, which is in better 
agreement with the criterion derived bv I Power et al. I l|2003l ). 
However, the results for these runs do not differ significantly 
from those with smaller softening length, so they are not 
presented here. 



2.3 Cosmological initial conditions 

From a cosmological 64^ h~^Mpc box we resimulate a halo 
promising to host a Milky Way-type galaxy. The crucial se- 
lection features are the halo mass (Mhaio — lO^'^ h~^ Mq), a 
relatively quiet merging history (the redshift of the last ma- 
jor merger with mass ratios of 1:3 and 1:10 are z=3.63 and 
z=2 respectively), a high spin parameter (A — 0.057) and no 
halo with equal or larger mass within a distance of 2 Mpc. 
We perform runs with three levels of resolution in the res- 
imulated region corresponding to an effective resolution of 
512^, 1024^ and 2048^. This corresponds to particle masses 
for dark matter (gas) of 1.47x 10** Mq (S.IxIO^Mq), 1.86 x 
10'' Mq (3.71 X 10*^ Mq) and 2.32 x 10*^ Mq (4.64 x 10^ Mq) 
respectively. Star particles have half the mass of the initial 



It is our goal to isolate the root cause of and find a solution ^ A halo with the density profi le as in equation [TO] quantified by 
for the angular momentum problem. In a fully cosmological iNavarro. Frenk fc White! l ll997^ . 
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gas particle mass. A gravitational softening of 4 (3), 2 (1.5) 
and 1 (0.75) li~^ kpc is used. 



2.4 Analysis of a simulated galaxy 

We define our halos using the virial radius Rvir, which is 
the radius where p(r) — Apc with pc = SHq/{8'kG) and 
A = 187r2+82a;-39a;^ x = Qo(l + zf /(no(l + zf + D.a)-1. 
The cold gas (T < 3 x 10*K) in the center of the halo makes 
up a clearly defined disk. In runs without star formation and 
feedback, all cold gas ends up in that disk. In the runs with 
feedback, some hot gas surrounding the cold gas is present, 
but the disk can still be distinguished clearly. When star for- 
mation is included, stars within the disk radius are defined 
as being part of the galaxy. Since the rotationally supported 
gas disk gives a clear axis of symmetry, we rotate the system 
to have this axis aligned with the axis in z direction and the 
disk in the x-y plane. The total specific angular momentum 
j of the disk is then computed by summing over all disk par- 
ticles, with respect to the center of mass of these particles. 
It is given by 



J = 



where 



7 



■mt{~T'i ~ ^com) X {'v'i 



V COM ) 



(7) 



(8) 



To compute the time evolution of the angular momentum, 
all particles belonging to the galaxy at z—Q are tracked back 
in time and the angular momentum is calculated in earlier 
snapshots for these particles with respect to the center of 
mass at the time. When star formation is included, the pro- 
genitor gas particles are used at the time when stars have 
not formed yet. 



3 THE ANGULAR MOMENTUM PROBLEM 
IN SIMULATIONS WITH RADIATIVE 
COOLING 

Our first, most basic set of runs includes dark matter and 
gas and the latter is subject to radiative cooling only. Star 
formation and feedback are turned off. This should enable us 
to pinpoint the processes responsible for the loss of angular 
momentum and investigate the role resolution can play to 
avoid this process without the additional, probably positive 
effects feedback should have. 



3.1 Isolated initial conditions 

We first analyze the results of our controlled semi- 
cosmological, isolated experiment in four different resolu- 
tions. Some characteristic parameters of the different runs 
are summarized in Table [T] 



3.1.1 The resulting disk 

The evolution in the semi-cosmological setup is illustrated 
in Figure [T] for the 32^ resolution. Since it is a physical co- 
ordinate system, initially the sphere expands following the 
expansion of the universe. The redshifts of turnaround at 




Figure 1. The formation of the disk in our 32^ sphere from the 
starting redshift of 75 until z=0. Each plot has a side length of 
280 physical kpc except the last one which zooms in to 80 kpc 
side length. 



z~12 and collaps e at z~8 agree wel l with analytical expec- 
tations based on iLokas fc HoffmanI l|200ll) . Small, cold gas 
clumps merge to form a central disk. This process is finished 
around z=l after which there is no longer any significant 
evolution. There is no late infall owing to the missing large 
scale environment. The disk can be described as a central 
overdense clump with a radius of about 5 kpc containing up 
to 80% of the total disk mass, and a thin outer disk. The 
resulting disk is unchanged with increasing resolution. 

We examine the s tability of the disk with respect to the 
Toomre Q parameter l|Toomrelll963 ) 



Qir) 



(9) 



where c^ is the effective soundspeed, k the epicyclic fre- 
quency and Eg the gas surface density. The disk is stable 
where Q(r)>1.0. As illustrated in Figure (2] in our case the 
disk is clearly unstable in the innermost, compact region, 
and only marginally stable in the outer parts (solid curve). 
Star formation can therefore happen in the disk. It could 
be stabilized if the gas were to remain hot with a minimum 
temperature of 5 x lO** K, as shown by the dashed line in 
Figure [1 

3.1.2 Verification of the setup with an analytical model 

We attempt a verification of our semi-cosmological initial 
conditions setup by comparing the resulting disk to predic- 
tio ns obtained with th e analy tical galaxy formation model 
bv IMo. Mao fc Whitel (|l998l ) (hereafter MMW). In their 
model MMW derive the properties of disks (like mass, scale 
radius, central surface density) from the density profile of 
the dark matter halo and its angular momentum under the 
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16^ 


32^ 


64^ 


128^ 


cosmological run 


DM particle mass [Mq] 


3.62x10** 


4.5x10^ 


5.64x10*5 


7.05x10^ 


1.86x10'^ 


gas particle mass [Mq] 


7.24x10'^ 


9.01x10'^ 


1.13x10*5 


1.41 xlO^ 


3.71 xlO** 


# of gas particles in disk 


1510 


12710 


106319 


606895 (z=1.7) 


38856 


e[h-l kpc] 


1 


1 


1 


0.25 


1.5 


Rvir [kpc] 


248 


251 


251 




292 


Vvir [km s-1] 


123.38 


124.45 


124.55 




145.45 


Mvir [Mq] 


8.73x10" 


8.99x10" 


9.01x10" 




1.21x10^2 


C 


28.04 


29.68 


30.542 




10.13 


A 


0.083 


0.088 


0.083 




0.055 



Table 1. Characteristic parameters for the semi-cosmological simulations with increasing resolution and for the cosmological run with 
an effective resolution of 1024'^ particles in the high resolution region as used in the run without star formation and feedback. 




Figure 2. Disk stability is measured with the Toomre Q parame- 
ter in our 32'^ sphere. The solid horizontal line shows the stability 
criterion Q > 1. The solid curve is for the standard model, the 
dashed curve for the model with 5 x 10* K minimum gas temper- 
ature. 



following four assumptions: (1) the disk mass is a fixed 
fraction of the halo mass (md=Mdisk/Mhaio), (2) the disk 
angular momentum is a fixed fraction of the halo angular 
momentum (jd=Jdisk/Jhaio), (3) the disk is thin, centrifu- 
gally supported and has an exponential surface density pro- 
file, (4) only dynamically stable systems can correspond to 
real galaxies. The disk properties are then completely de- 
termined by the rotational velocity of the halo at r2oo (the 
radius where p{r) = 200pc), V200, the concentration c of the 
halo, the spin parameter A, md and jd, and are calculated 
iteratively taking into account the influence of the forming 
disk on the halo density profile and rotation curve. This is a 
gravitational infiuence as well as a contraction of the inner 
regions of the halo. 

As input parameters we use the results from our simulations, 
the disk mass Md and disk angular momentum Jdisk from 



400 



300 - 



200 




20 30 
R [kpc] 

Figure 3. Comparison of the rotation curves from the simulation 
(red) and predicted by the MMW model (black). 



the run with radiative cooling, and the dark matter halo 
parameters (halo mass M200, halo angular momentum J200, 
halo concentration c, spin parameter A) from the pure dark 
matter run (note the use of R200 instead of Rvir here, to 
define the halo in agreement with MMW) . For our purposes 
we mainly look at the predicted rotation curves for the disk 
and dark matter halo, which we then compare to the rota- 
tion curves obtained in our simulations. Table [2] gives the 
MMW parameters for the different runs. We applied two 
major modifications to the original MMW model. MMW 
assume an NFW halo with a density profile given by 

Pir) = Pc , / w, 1 , S2 ' (10) 
(r/rs)[l + r/rsY 

where 5c is a characteristic density and is a scale length. 
Since our initial conditions are only semi-cosmological and 
we are missing the surrounding large scale structure, we are 
also missing some late infall of material into our halo which 
would happen in a fully cosmological box. This results in 
density profiles which deviate from a typical NFW profile at 
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large radii. Therefore, instead of using a halo mass profile 
based on the NFW profile as done by MMW, we use the 
calculated mass profile of the dark matter halo as an input 
for the model. 

A second problem we encountered was a systematic 
over-prediction of the inner part of the dark matter rotation 
curve by the MMW model. As already mentioned above, 
the mass profile of the halo is influenced by the formation 
of the disk. This effect is describe d by the adiabat i c con - 
trac tion model, based on wo rk by I Barnes fc White! (|l984l ) 
and iBlumenthal et al] (Il986l ). In this model, particles are 
assumed to move on circular orbits and angular momen- 
tum is strictly conserved. For circular orbits, the force is 
F/m = n^/r and the angular momentum is j /m = vr. Com- 
bining them and including angular momentum conservation 
leads to the general result for a particle starting at radius ri 
and ending at radius rf 



f 



Fin, 



(11) 



With the Newtonian force law F/m= ^'^{j^ ) ^ this leads to 
Mf{rf)rf = Mi{ri)ri, where M(r) is the enclosed mass at 
radius r. This is the form used by MMW. However, in or- 
der to reduce the occurrence of two-body direct interactions 
in the simulations, the forces are softened and may be de- 
scribed rather with a Plummer force law 
_ ~GMmr 

Fpiummer = ^2)3/2 ' (l^) 

We therefore modify the MMW model by using the force 
profile inste ad of the mass pr ofile to calculate the adiabatic 
contraction (|Gottbrathll200lh . The force profile is given by 



(13) 



Fir) = / F{S,r)dS 



where F(S,r) is the Plummer force profile of a spherical shell 
with radius S which in turn is given by 



dF _ GmM / + Sjr + S) _ + SjS - r) 
dS ~ 2r^S \ + [r + 5)2 + [r - S)^ 



.(14) 



With this modified approach we achieve good agreement 
between the rotation curves predicted from the model and 
produced in the simulation, as can be seen in Figure [S] The 
only region being different is the inner 15 kpc of the disk ro- 
tation curve, which is much more peaked in the simulations. 
This is due to the high concentration of matter in the center 
as a result of the angular momentum problem, and to a lack 
of resolution in the center. The central clump behaves like 
a bar, sweeping up gas directly surrounding it and resulting 
in an underdense shell between the central clump and the 
thin disk. This is also visible in the bottom right plot of the 
disk in Figure [1] Therefore, the rotation curve is too peaked 
in the very center, and too low between radii of 5 and 18 
kpc. 



3.1.3 The angular momentum analysis 

We analyze the loss of angular momentum by calculating in 
all snapshots the specific angular momentum of cold gas par- 
ticles within a radius of 250 kpc (which is the virial radius 
at z=0). Since most cold gas particles end up in the disk, 
this basically traces the angular momentum evolution of disk 
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128" 
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Figure 4. The evolution of the specific angular momentum over 
time for the cold gas within the virial radius. The different colored 
lines show the 4 different resolutions of the semi-cosmological ha- 
los. The dash-dotted and dash-double-dotted parts of the blue 
line indicate the evolution for the hightes resolution case, when 
gas was converted to stars, calculated in two different ways. As a 
comparison, this was also done for the 64^ resolution, shown in 
the green dash-double-dotted line. 



particles. The result is shown in Figure |3] The highest res- 
olution run (128^ particles in the underlying box, blue line) 
poses a computational challenge. Without star formation, it 
is computationally infeasible to advance the simulation be- 
yond z~1.7. We overcome this by replacing all gas particles 
at this point with coUisionless "star" particles. In this way, 
no further cooling is possible, but the disk is already largely 
in place at this point. Reducing the computation to only 
gravitational forces enables us to extend the run until z=0. 
For the angular momentum calculation unilt z=1.7, we use 
the cold gas within the virial radius, as in the other resolu- 
tions. After that we follow two different methods. First we 
tested to use all star particles within the virial radius (thin, 
dash-dotted blue line). However, in a comparison of the same 
procedure for the 64^ resolution, we find much better agree- 
ment with the original gas run when following only the stars 
making up the galaxy at z=0 (all star particles in a radius 
of 50 kpc). This is shown in the thin green dash-dotted line. 
Applying the same procedure to the 128^ resolution results 
in the lower, blue dash-triple-dotted line. The discontinuity 
to the gas run at z=1.7 is due to the switch in calculation 
methods. The resulting angular momentum evolution of the 
two calculation methods are not very different. The result 
of the real gas disk is expected to be in between the two. 
In both cases, there is good agreement with the evolution 
in the lower resolution runs, strongly indicating a gravita- 
tional effect as the reason for the angular momentum loss 
and ruling out artefacts from the hydrodynamics or SPH 
implementation in the code. The behavior in all four reso- 
lutions is quite similar. For the highest resolution the loss is 
slightly stronger, indicating that resolution makes the prob- 
lem more severe instead of improving it. We first calculate 



The Angular Momentum Problem in Cosmological Simulations of Disk Galaxy Formation 7 



resolution 


R200 


V200 


M200 


j200 


disk 


jdisk 




jd 




[kpc] 


[km s^-^] 


[M©] 


[kpc km s^-"-] 


[Mq] 


[kpc km s^-"-] 






163 


189 


138.24 


8.36x10" 


2233.25 


1.09x10" 


1017.6 


0.13 


0.059 




191 


139.42 


8.59x10" 


2100.72 


1.15x10" 


1092.4 


0.134 


0.07 


643 


192 


139.47 


8.64x10" 


2230.24 


1.2x10" 


1234.8 


0.139 


0.077 



Table 2. Parameters for the MMW model. 



the ratio of the final over the initial specific angular momen- 
tum of the disk particles for our four resolutions. Over our 
entire resolution range we find the same net loss of angu- 
lar momentum of about 55% and no trend with resolution 
is observed. This indicates that a simple increase in resolu- 
tion does not help to alleviate the angular momentum prob- 
lem within the framework of the basic hierarchical model of 
galaxy formation. 

We now investigate, when and how the angular momen- 
tum is lost. In combination with Figure [T] Figure |4] clearly 
illustrates that the bulk of angular momentum is lost during 
the phase of evolution dominated by the infall and merging 
of smaller clumps. At z=l, the angular momentum ratio 
is ja^i/jinit ^ 0.55. Once the disk is established at z~l, 
there is still some loss, but this happens gradually with 
jfinai/jz=i ~ 0.88. This further illustrates that improved res- 
olution does not strongly influence this process, since the 
loss not only follows the same pattern in all resolutions, it 
also would otherwise continue stronger even after the forma- 
tion of the disk. We can further str engthen our conclusion 
by foUowing lWeil. Eke fc Efstathioul ((1998 ) and evolving the 
gas component adiabatically until z=l. The lack of radia- 
tive cooling prevents the formation and therefore merging 
of small clumps in the gas. The disk forms quickly after the 
coohng is turned on at z=l. Figure [5] shows the resulting 
difference for the gas in the 32^ run. There is still a loss of 
angular momentum, but it is smaller than before, with a final 
ratio for the disk of jflnai/jinit=68.3% and jz=i/jinit=95.1% 
at z=l. 



3.2 Cosmological initial conditions 

We now turn our attention to the angular momentum prob- 
lem in a fully cosmological run to verify our conclusions. For 
this we use the medium resolution cosmological halo, with 
1024^ particles effective resolution in the high resolution re- 
gion. Characteristic parameters of the halo and disk can be 
found in Table [T] 

The disk in this case is as centrally concentrated as in 
the semi-cosmological case with 78% of the disk gas within 
the central 5 kpc. The build-up of the gas disk is shown 
in Figure |S] The top and middle panels have a side length 
of 1 h-i Mpc in comoving coordinates and nicely show the 
hierarchical growth from small clumps to larger objects. In 
the time between the two middle panels (z~1.5 and z~0.75) 
a last larger merger of two gas clumps occurs. After that, 
the disk grows only via accretion. The bottom panel zooms 
into the disk before and after that merger. It is interesting 
to now compare (see Figure [7]) the mass growth history of 
the disk (solid line) and the evolution of the specific angular 
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Figure 5. Comparing the angular momentum evolution for the 
32^ halo in the standard case (with radiative cooling switched 
on all the time, black thick curves) with the run with adiabatic 
evolution until z=l and cooling after z=l (blue thin curves). The 
solid lines are for disk gas, the dashed lines for all gas within the 
virial radius at z=0. 



momentum of the disk particles (dashed line). The arrows 
indicate the points in time at which the snapshots in Fig- 
ure |6] are shown. The early growth of angular momentum 
accompanying the first collapse of the halo follows well the 
expectations from tidal torque theory (Joct). The maximum 
angular momentum acquired by the gas agrees well with 
the angular momentum of the dark matter halo. However, 
it can clearly be seen, that in particular the merger of two 
large gas clumps at z«1.5, which increases the mass of the 
disk by about 25%, is accompanied by a large loss in angu- 
lar momentum (about 50%). The ratio of final to maximum 
specific angular momentum is 0.37. This is a little less than 
what we found in the semi-cosmological case, as expected, 
thus confirming and even strengthening our conclusions of 
the previous section. 



4 THE INFLUENCE OF STAR FORMATION 
AND STELLAR FEEDBACK ON THE 
ANGULAR MOMENTUM PROBLEM 

In Figure [8l we show the gas and stellar component of the 
resulting galaxy in our simulations including star formation 
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Figure 8. Results at z=0 of the runs with star formation. In the first and third columns we show runs without feedback, and runs with 
feedback are displayed in the second and fourth column. The left two columns are semi-cosmological results, the right two columns are 
fully cosmological runs. The face on and edge on plots have a side length of 80 kpc. The top two rows show the cold gaseous disk, rows 3 
and 4 the stellar component. The star particles are color-coded by age (red: older than 10 Gyr, blue: younger than 4 Gyr). The bottom 
row shows the star formation history. 
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Figure 7. Growth of the disk (solid hne) and evolution of angu- 
lar momentum (dashed line) for the fully cosmological run. The 
big loss in angular momentum starting at z=1.5 clearly coincides 
with a large mass increase. The arrows indicate the time of the 
snapshots shown in Figure |6] 



noticeable, that the gas disk is much more extended than the 
young stellar component in the cosmologica l runs. This is in 
agree ment with observations (for example iBoomsma et al.l 
I2OO8II . 



Figure 6. Snapshots illustrating the assembly of the gas disk in 
the fully cosmological runs. The top four plots have a side length 
of 1000 comoving kpc h"'^ to show the larger environment, the 
bottom two plots zoom in to show more clearly the disk itself and 
has a side length of 200 kpc h~^ in physical coordinates. 



(columns 1 and 3 for the semi- and the cosmological en- 
vironments, respectively) and including star formation and 
feedback (columns 2 and 4 for the semi- and the cosmological 
environments, respectively) . The first feature to notice is the 
ring-like structure of the gas disk in the semi-cosmological 
runs (top row, columns 1 and 2). The reason for this is the 
efficient transformation of gas into stars moving inside-out, 
and the missing late infall of gas to replenish the disk. The 
disk here is also much more defined than in the fully cos- 
mological case, again due to the missing infiuence of late in- 
coming clumps. In all cases, the young stars (younger than 4 
Gyr, colored in blue) and the cold gas (T < 3 x 10^ K) form 
a fairly thin disk. In both the semi-cosmological and the 
fully cosmological run without feedback (column 1 and 3) 
stars are formed very efficiently at early times, resulting in a 
dominant spheroidal component. This can also be seen in the 
strong peak in the star formation history at a t~2 Gyr (bot- 
tom row). When feedback is included (column 2 and 4), star 
formation is quenched efficiently and a more prominent disk 
of young stars is formed while the large spheroidal compo- 
nent made up of old stars (colored red) is reduced. It is also 



4.1 The origin of the stellar components 

Why do we find such a large spheroidal component in the 
runs without feedback, even though the gas disk in the 
runs without star formation is very thin? The answer can 
be found in the location and the time where the stars are 
formed. In our fully cosmological run, star formation starts 
at z~9, and it peaks at z~3. This is the time when the 
progenitor of the z=0 disk is just starting to be assembled, 
so the majority of the stars are actually forming in larger 
mass clumps which fall in until z~1.2. After that, mass 
grows more slowly through accretion and the infall of smaller 
clumps. Stars are continuously formed in the disk progeni- 
tor, but at a very low rate of about 2 M0 yr"'^. The stars 
entering the halo in clumps are not able to settle into a disk 
and therefore will form a spheroid. Only the stars formed 
within the disk progenitor will be rotationally supported 
and can form a stellar disk. This can be seen in Figure |9l 
which shows the distribution of stars at z=0, divided accord- 
ing to their formation environment into "clump stars" (stars 
formed outside of the disk progenitors and assembled into 
the galaxy via merging) and "disk stars" (stars formed in the 
progenitor of the z=0 disk). "Disk stars" form a reasonably 
thin distribution tracing the gas disk, "clump stars" make 
up the large spheroid. One important diagnostic of the na- 
ture of the stars is the ratio of rotational velocity to velocity 
dispersion. Both parameters are plotted in Figure [10] for all 
stars (black line), "clump stars" (red line) and "disk stars" 
(blue line). For "disk stars", this ratio is much higher, con- 
firming their rotational support compared to dissipational 
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Figure 10. Rotational velocity (top) and velocity dispersion 
(bottom) of all stars (black, solid line), "clump stars" (red, dashed 
line) and "disk stars" (blue, dash-triple-dotted line). 



Figure 9. Face-on and edge-on projection of the stars at redshift 
0. The left column shows the stars which were formed in clumps 
and then accreted. The right column shows the stars which were 
formed in the disk progenitor. 



support for "clump stars". Of all stars in the galaxy at z=0, 
30% are "disk stars" and 70% are "clump stars". We per- 
formed a test where star formation was only turned on at 
z=l, with the thin gas disk being more or less in place and 
further mass infall mostly through smooth accretion. The 
stars in this case indeed form a very prominent stellar disk, 
but no large spheroidal component. 

In the feedback runs, star formation also begins early 
in individual clumps. The feedback however heats the gas 
and pushes it out of these structures, quenching the star 
formation in the clumps and enabling the gas do be ac- 
creted smoothly later directly into the disk where it then 
can form disk stars. This is evident from the reduced early 
star formation peak in the feedback runs as shown in Fig- 
ure |8] Consequently, in the final galaxy we have only 35% 
of all stars formed in clumps (and therefore a much less 
prominent spheroidal component as visible in Figure |8] right 
column) and 64% formed in the disk progenitor. 



ing the cold gas and stars of the final galaxy at z=0 over 
time. In all three cases we see a drop in angular momentum 
at z~1.2, but it is less pronounced in runs with star forma- 
tion and even smaller in those with feedback. However, in 
simulations with star formation there is some loss of angular 
momentum after z=l, which is not as significant when only 
cooling is employed. This is probably due to the increased 
late accretion of gas and related hydrodynamical torques. 
Compared to the maximum angular momentum, an even 
slightly higher percentage is lost in the case with star for- 
mation, and still about 50% with feedback. Most of this loss 
occurs after z=l. The feedback has reduced the loss owing 
to dynamical friction but increased the later loss due to hy- 
drodynamical torques. While a loss of 50% still seems large, 
this is consistent with exp ectations comparing dark m atter 
halos and observed disks (|Navarro fc Steinmet j|200oh . For 
the feedback case, we plot the results of runs with three dif- 
ferent resolutions (effective resolutions of 512'^, 1024^ and 
2048''), which behave similarly indicating that resolution 
does not affect this process much. We conclude that feedback 
in both the semi-cosmological and the cosmological simula- 
tions clearly helps to alleviate the angular momentum prob- 
lem which is mainly due to the reduction of the early star 
formation and the resulting smaller spheroidal component. 



4.2 The angular momentum problem 



5 SUMMARY AND CONCLUSIONS 



Since now gas particles are converted into stars, we compute 
the angular momentum for the combination of cold gas and 
stars. The result for the semi-cosmological runs is shown in 
Figure |ll(a)[ The loss is smaller when including star for- 
mation (dashed line), but still substantial. In the case with 
feedback (dotted-dashed line), later accretion of gas with 
high angular momentum keeps the angular momentum of 
the disk high and the final galaxy has almost 90% of the ini- 
tial angular momentum. Gas is available for late accretion 
due to the decreased early star formation. The results for the 



cosmological runs are shown in Figure 11(b) We are track- 



In this paper we have studied the causes of the angular 
momentum problem with special emphasis on the influ- 
ence of numerical resolution and the impact of star forma- 
tion and feedback. We performed a systematic study lead- 
ing us from a controlled, isolated, "semi-cosmological" ex- 
periment to fully cosmological simulations. We also inves- 
tigated simulations with different physical processes, once 
only taking into account radiative cooling, then including 
simple star formation and finally also including supernova 
feedback. The advantage of the semi-cosmological environ- 
ment is to combine at moderate computational costs an 
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Figure 11. These plots compare the specific angular momentum evolution in the runs without star formation (solid line), with star 
formation but without feedback (dashed line) and with star formation and feedback (dashed-dotted line). Panel (a) shows the semi- 
cosmological, panel (b) the fully cosmological runs. The three different lines for feedback in Panel (b) show runs with three different 
resolutions. 



isolated, well controlled halo with the small scale hierar- 
chical buildup which is characteristic for a ACDM uni- 
verse. By performing first simulations without star forma- 
tion, it is possible to assess the importance of resolution 
in the hierarchical framework for structure formation, es- 
pecially with respect to the angular momentum problem. 
Our results clearly show that an increase in resolution in 
this case does not help against the angular momentum loss. 
The loss of angular momentum for the disk gas is of the 
same order while the resolution changes by orders of mag- 
nitude. We are also able to tie the angular momentum loss 
to the phase of collapse into small cold gas clumps which 
merge and form the main object. The main loss mecha- 
nism is dynamic al friction of the clump s and larger merg- 
ers, as shown bvlNavarro fc Ben3 ( 199lf). iNavarro fc White! 
lll994l.lNavarro fc Steiiimetj (llOgTi v ld'Onghia et al.l (120061 ) 
and Izavala et al.l (|2008| ). The efficiency of dynamical fric- 
tion is dominated by the largest clumps. These c lumps can 
be tr acked even in low resolution simulations (|Gao et al.l 
I2OO4I . and refer ences therein). This ex plains why, contrary 
to for example iKaufmann et al] (|2007l ) , we find no signifi- 
cant resolution dependence even t hough our resolution range 
is comparable. Kauf mann et al.] (|2007l ) studied the forma- 
tion of a gas disk starting in an equilibrium halo without 
taking the halo formation history into account, so their an- 
gular momentum loss is mainly due to hydrodynamical and 
gravitational torques during gas infall. This is a resolution- 
dependent process. We also see this effect, as can be seen 
in the run with cooling switched on only after z—1 (Fig- 
ure [Sjl and in the cosmological runs with star formation and 
feedback. However, the effect is much smaller than the loss 
during the merger period. We therefore argue that, in order 
to prevent the angular momentum loss in cosmological sim- 
ulations, mainly the early formation of gas and star clumps 
must be prevented. The loss owing to limited numerical res- 
olution appears to be only of second order. In line with this 
argument, we observe a comparable loss in our fully cos- 
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Figure 12. For the semi-cosmological and cosmological run, we 
compare the final angular momentum of the disk in the cooling 
run ("cool"), with star formation but without feedback ("w/o") 
and with feedback ( "with" ) . The disk angular momentum is com- 
puted once including the cold gas disk and all stars in the galaxy 
(crosses), and also including the cold gas and the young stars only 
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mological simulation. It can be clearly linked again to the 
merging of two large gas clumps. Unfortunately it is compu- 
tationally not feasible to run this at higher resolution with- 
out star formation. However, even if resolution would have 
some infiuence in cosmological simulations (contrary to the 
semi-cosmological case), this would only make the problem 
worse by increasing the number of resolved gas clumps. We 
are confident that this would not lead us to a different con- 
clusion. 

The most plausible mechanism to prevent clump for- 
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mation is star formation and related supernova feedback 
because it heats the gas. Another possible source of heat- 
ing would be the cosmic UV background. Some low resolu- 
tion tests we performed inclu ding this effect did no t show 
any positive impact. Indeed, iNavarro fc Steinmet3 (Il997l ) 
showed that a UV background may even have an adverse 
effect as it suppresses the infall of gas which is accreted late 
in the formation history and thus acquired high angular mo- 
mentum. Here we show the changes in the angular momen- 
tum evolution when star formation and supernova feedback 
is included. As expected, star formation alone does not re- 
sult in an improvement, since the same clumps are spiraling 
in. They simply consist of highly condensed stars instead of 
highly condensed gas. Star formation peaks high early on, 
and especially in the semi-cosmological setup, where late gas 
infall is missing, most gas is used up at z—Q. Feedback helps 
to prevent this by keeping the gas too hot to form stars and 
therefore allowing it to settle into the disk later on as high 
angular momentum gas. We studied the origin of the stars 
in the final galaxy and found that without feedback 70% 
of all stars in the final galaxy enter in clumps. With feed- 
back this is reduced to 35%. This can also be seen in the 
angular momentum of the final galaxy. Clump stars form a 
slowly rotating spheroidal component and disk stars form a 
rotationally supported, high angular momentum disk. The 
final angular momentum of the gas and all stars combined 
is therefore much higher in our runs with feedback (and a 
smaller spheroid) than in the ones without feedback (and a 
large spheroid). This is further illustrated when looking at 
the combined angular momentum of the rotationally sup- 
ported component of the galaxy, consisting of the cold gas 
and the young stars only (Figure [12}. This is roughly equal 
in the runs with and without feedback (triangles), while for 
the total galactic angular momentum including all stars a 
clear increase can be seen from the run including only ra- 
diative cooling over the run with star formation to the run 
with star formation and feedback. This increase is due to 
the decrease of the central gas clump and then of the stellar 
spheroid. 

We conclude that an efficient suppression of the formation of 
dense clumps at early times is the main factor in reducing 
the angular momentum loss. Star formation alone cannot 
provide this effect since these dense gas clumps also fulfill 
the conditions for star formation resulting in a high early 
star formation peak and the formation of a large bulge. Only 
feedback can reduce this effect efficiently. Resolution over- 
all plays a secondary role. In our simulations with star for- 
mation and feedback, the angular momentum loss is about 
50% which is in agreement with expectations from compar- 
ison s between dark matter hal os and observed disk galax- 
ies (|Navarro fc Steinmet j[2"o00l ). This loss mainly happens 
at late times due to hydrodynamical torques. In our sim- 
ulations with three levels of resolution in the cosmological 
halo, resolution does not influence this effect either. While 
feedback has such a critical role in the formation of more re- 
alistic disks in simula t ions, it is still not well understood. In 
iPiontek fc Steinmet j (l2009bl ) we therefore present a detailed 
study of stellar feedback in disk formation simulations. 
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